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f^ Abstract. An analysis is provided of the degradation that arises in the quality 

£vj factor of a whispering gallery mode when a circular or spherical dielectric cavity is 

deformed. The large quality factors of such resonators are important to their use 

— i in applications such as sensors, wavelength filters or lasers. Yet a straightforward 

^_^ analysis of the effect of shape deformation on quality factors cannot given because 

^^ the underlying complex ray data demanded by a standard eikonal approximation 

frequently does not exist. In this paper we exploit an approach that has been 
C*l successfully used elsewhere to describe the strong directional emission of such 

CN systems, based on a perturbative treatment of the relevant complex ray families. 

Applicable when the radial perturbation is formally of the order of a wavelength, 
the resulting approximation successfully describes changes to the quality factor 
using the ray geometry in a neighbourhood of a discrete set of escaping rays 
guiding the directions of maximum emitted intensity. 
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Optical microcavities are widely exploited in applications such as lasers, sensors and 

wavelength filters [1, 2, 3, 4]. A key underlying feature in such technologies is that 

<**> | these devices support whispering gallery modes, which, corresponding to ray families 

confined by total internal reflection to the interior of the resonator, are very long lived. 

,__! An important practical consideration is then to be able to characterise the quality 

^ factor of such systems. Geometry is a particularly important factor in this regard. 

Deformation from perfectly circular or spherical geometry, while leading to the often 
desirable feature of directional emission, typically degrades the quality factor itself 
[5, 6]. Here we treat the case of smooth shape deformations but note that qualitatively 
similar features may also be seen in microcavities perturbed by notches [7] or internal 
scatterers [8, 9, 10]. 

In this paper we focus particularly on the regime of weak deformation and offer 

^-H analytical models for the quality factor in that limit. Whether deliberately engineered 

or arising as a consequence of manufacturing tolerances, even very weak deformations 

. — can be be seen to have a dramatic effect on the external field of these systems and as 

a consequence (we will find) on the quality factor. Despite the fact that the internal 
modes and associated real ray dynamics are hardly qualitatively changed from the 
corresponding exact circular or spherical limit, a straightforward application of eikonal 
theory frequently fails for such systems because natural boundaries [11, 12, 13] may 
prevent one from calculating the complex ray data required by an approximation of the 
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external, evanescent field. Here we exploit an approach which has successfully been 
used elsewhere to characterise the directional emission of such systems [14, 15]. Rather 
than demand exact solutions of the eikonal equation governing ray dynamics, which 
natural boundaries may prevent us from continuing sufficiently far into the complex 
domain, we instead construct approximate ray families, calculated using canonical 
perturbation theory. 




Figure 1. The external field is shown for a whispering gallery mode in an elliptical 
cavity with refractive index n = 2 and focal half-length c = ro/4 where ab = r^ is 
the area. This mode has azimuthal quantum number m = 50 and 5 radial nodes 
(with k r TQ fa 37.5) and we note that the external field is exaggerated graphically 
to make it visible. Superimposed are the real escaping rays 71 and 72 guiding 
waves of maximum intensity and the real rays 73 and 74 guiding waves of minimum 
intensity. The geometry of the ray family around these particular escaping rays 
determines the Q factor. 

Once the external field has been successfully approximated in this manner, the 
quality factor itself may be calculated using Herring's formula [16], which represents 
it as a flux integral measuring the overall rate of escape from the resonant mode. 
A straightforward evaluation of this flux integral, following a truncation to first 
order of a perturbative approximation of ray families, is shown in explicit numerical 
examples to describe significant degradations of the quality factor in problems where 
the deformation seems at first sight to be very weak. It should be emphasised 
that, although the underlying ray calculation is perturbative, it successfully described 
emitted waves that are strongly altered by the deformation and applies when 
perturbation of the wave solutions itself is not valid [17]. 

When approximated for somewhat larger deformations using saddle-point 
approximation, one finds that the dominant contribution to the Herring integral is 
associated with a discrete set of rays escaping to infinity along the directions of greatest 
emitted intensity (such as 71 and 72 in Fig. 1). Dynamically, the distinguishing feature 
of these rays is that they propagate with real direction cosines, whereas one finds 
generically in deformed resonators that escaping rays are slightly complex. In terms 
of phase-space geometry, the escaping ray family locally forms a complex Lagrangian 
manifold and the dominant real escaping rays form the intersection of this Lagrangian 
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manifold with its complex conjugate. For sufficiently large deformations, the leading 
order approximation of the quality factor is determined by a variant of Wilkinson's 
formula [18, 16], which expresses it in terms of canonically invariant measures if the 
transversality of this intersection. 

Although a primitive application of the steepest descents approximation fails in 
the limit of very small perturbation, we also offer a uniform analysis allowing one to 
interpolate smoothly from the case of small deformation, where simpler perturbative 
results apply, to larger deformations where Wilkinson's formula can be used. This is 
shown in particular to give a good description of the quality factor of elliptical cavities, 
where the existence of exact invariants allows us to evaluate the relevant complex ray 
data for a wide range of deformations. 

We conclude by summarising the content of the paper. We begin in Sec. 2 
by setting out the basic problem to be solved and the notation used. We also 
describe in general terms how Herring's formula can be used to calculate the quality 
factor and outline the main approaches to its approximation in the cases of weak 
deformation, where completely perturbative results are appropriate, and of moderate 
deformation where the Herring flux integral may be tackled using approximation by 
the method of steepest descents. A more detailed description of these approaches, 
and their results, are offered in subsequent sections. Sec. 3 describes the perturbative 
approximation which may be applied to weak but generic deformations. Sec. 4 
describes the asymptotic development of the Herring integral appropriate to larger 
deformations and applies this approach in particular to the special case of elliptic 
cavities, where the requisite complex ray data may be explicitly characterised, even 
for larger deformations. Finally, conclusions are offered in Sec. 5. 

2. Notation, Herring's formula and its approximation 

2.1. Notation 

We consider a two-dimensional scalar model governed by the Helmholtz equation 

-V 2 -</> = n 2 (x)fcV- 

(We expect that the basic underlying approach offered here, and the general 
conclusions reached, will apply also to higher dimensional problems or vector 
formulations, but we restrict our detailed discussion to two-dimensional scalar waves 
for clarity of exposition.) Denote the region occupied by the cavity by fl and its 
boundary by dfl, so that the refractive index changes discontinuously from a constant 
value n>linf2ton=l outside O. The resonances we treat are examples of 
Gamow-Siegert states [19, 20], and are solutions of the Helmholtz equation subject 
to outwards-radiating boundary conditions at infinity. Such resonances are found for 
discrete complex wavenumbers 

2 
with a negative imaginary part, which is small when the resonances are long-lived, as 
is the case for the whispering-gallery modes we consider. The quality factor 
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is then large. The Gamow-Siegert states are not normalisable and instead we adopt a 
normalisation convention 

m 2 a = /i^( X )i 2 dx = i, a) 

in which integration is restricted to the cavity's interior. 

2.2. Herring's integral 

Herring's formula is obtained by integrating the identity 

over a region R containing CI, and using Green's identities to get 



2i£v||?/>||jj. j E \ dn dn 

where s represents arc- length along £ = dR and where HV'lllj is defined as in 
(1), but with integration restricted to R instead of CI. The result is, of course, 
independent of the choice of region R, although some restriction on this choice makes 
the approximation that follows easier. In particular, we choose the section £ to 
be far enough outside Cl that the rays underlying an eikonal approximation of ip 
are approximately real (see below), but close enough to Cl that the exponentially 
small contribution to \\tp\\% from the exterior of Cl can be neglected, allowing us to 
approximate 

The Herring integral (2) then measures the flux of the quasinormalised state ip escaping 
across the section £. We now outline how an asymptotic approximation of it is 
determined in terms of the geometry of the ray family underlying ip. 

2.3. Approximation of the wavefunction 

In the case of whispering gallery modes in a circular cavity, the wave field immediately 
outside the resonator decays exponentially in the radial direction and is approximated 
by a family of rays whose radial momentum is purely imaginary. Far enough outside 
the resonator, however, a caustic is encountered beyond which the rays once again 
propagate with real radial momentum and the wavefunction is oscillatory, albeit 
with an exponentially small amplitude following the decay that occurs through the 
evanescent band immediately outside the resonator. We consider the case of deformed 
resonators sufficiently close to circular that a deformation of this ray family can be 
defined which provides an eikonal approximation 

V^(x) rj A(x)c ifeS ( x > 

for the field in the exterior region. Although natural boundaries may prevent us from 
defining such a ray family exactly, we have shown elsewhere that approximate families 
can be constructed using canonical perturbation theory which successfully describe 
the deformed exterior wave field [14, 21]. We explain this construction in more detail 
later but for now we simply assume that such a ray family can be defined, with an 
overall geometry that is a deformation of the circular case. 

We consider in particular the case where x is far enough outside the resonator 
that the rays of the circular limit are real. In the deformed case, a symmetry breaking 
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is found to occur so that the corresponding rays become slightly complex. That 
is, even for x G R 2 , the rays propagate with complex direction cosines. Then the 
resulting direction dependence of the imaginary part Im(S'(x)) of the action function 
can strongly modulate the wave intensity |?/;(x)| 2 and one finds that the Herring 
integral may be dominated by a few directions of maximal flux density. We now argue 
that the condition of maximum intensity/flux corresponds in terms of ray geometry 
to the condition that the escaping rays are real, and provide a general description of 
the asymptotic evaluation of the flux integral around these real ray directions. 

2.4- Approximation of the Herring integral 

To describe the asymptotic evaluation of the Herring integral (2) more explicitly, let 
us write the amplitude of the wave field in the form 

where t is a transmission coefficient coupling the interior and exterior solutions, defined 
so that the remaining density term takes a standard WKB form. [We point out that 
a full WKB treatment of dielectric cavities requires the inclusion of a phase — related 
to the Goos-Hanchen shift [22, 23] — accounting for variation of the reflection phase 
along the boundary. To simplify the presentation as much as possible we assume 
here that any such phases have been incorporated into the transmission amplitude t.] 
Then, at leading order (in an expansion in powers of 1/fc), Herring's integral becomes 

ft ' (^^r-) W ds « / (RcPn)\t\ 2 \ P \ e- 2MmS ds, (3) 



'£ \ * / JY, 

where n is the unit normal to £, p(x) = VS 1 and p n = n • p. 

Further simplification of this integral can be achieved by exploiting the associated 
dynamical invariants for modes with an intcgrable underlying ray dynamics. Any ray- 
dynamical system enjoys the invariant H = (p 2 +p 2 )/2. An additional invariant, which 
we denote by M, is assumed to exist for the ray families underlying the modes we treat. 
This second invariant is angular momentum for the circular limit and a deformation of 
it for perturbed cavities. In the case of elliptical deformations, an explicit analytical 
form will be given for M in Sec. 4 whereas for the case of generic deformation, an 
approximate invariant M is defined implicitly by the canonical perturbation approach 
we employ in Sec. 3. Let \ denote an angle coordinate on £ conjugate to M. That is, 
X evolves at a constant rate under the flow generated by using M as a Hamiltonian 
and we scale it so that it has period 27r. Then, it can be shown that 

J dx 



\p\ ds ., I 

27T \p n \ 



where 

J 



d(<Pi,<P2) 



d(t, X ) 



d(H 7 M) 



d(Ji,J 2 ) 

is a Jacobian relating the variables (t, x), respectively conjugate to (H,M), to 
regular action-angle variables (cpi, (p 2l J\, J 2) for the ray family inside the cavity. The 
important feature for now is that J is constant on ray families. Herring's integral can 
now be written 

J f Rc P „ |i|2 e _ 2ferImS ^ (4) 



27Tj s \p n \ 
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Significant further approximation of this integral is possible in two particular regimes. 
First, in the case of small generic deformations on a scale e = 0(l/k), the 
imaginary part of the action can be computed perturbatively while the amplitude 
variation can be neglected and the quantities p n and t in this integral replaced by the 
constant values they take for the circular limit. Then 

K « K M, (5) 

where 

M = li c -2Wm(S( X )-S„) dX] (6) 

27T / 

where k$ and S'o(x) respectively denote the value of n and the action function for 
the circular limit. Further development of this approximation will be offered in Sec. 3 
where approximations for Im(S(x)) are developed by applying canonical perturbation 
theory to the ray families. We also note that the integrand e - 2k ^MS(x)~So) p r0 vides 
a leading-order measure of the intensity of the emitted field and has already been 
calculated in that context in [14]. The simple perturbative procedure pursued in this 
paper applies in the generic case where the unperturbed ray family is not too close 
to resonances of the ray dynamics. When ray-dynamical resonances are nearby, a 
somewhat more involved approximation of the exterior field is required but even in 
that case we note that the Q-factor degradation can be evaluated once the exterior 
field is known using 

M = J m 2 d X / J Ho\ 2 d x 

where ipo and ip £ respectively denote the circular limit and the ^-deformed form of 
the exterior field. Explicit calculations are confined in this paper to the simpler case 
described by (6), however. 

For larger deformations, it is natural to approximate the integral using the method 
of steepest descents. For generic deformations, however, the intervention of natural 
boundaries means that we often cannot find the required ray data. An exception is the 
case of elliptical cavities where a complete analytical description of the external ray 
family can be achieved. The resulting approximation for this case is described in detail 
in Sec. 4 but here we summarise the main qualitative features, which also apply to the 
large-efc asymptotics of (6). Stationary points of the exponent — 2Im(5) = i(S — S*) 
in Herring's integral are obtained for real x when 

P = P*, (7) 

that is, when the associated escaping rays are real. Generic escaping rays are complex 
for deformed resonators (albeit only slightly complex when the resonator is weakly 
eccentric). A discrete subset of real escaping rays can typically be found, however, 
propagating along the directions of greatest and least emitted intensity. A steepest- 
descents approximation of Herring's integral therefore characterises the decay rate in 
terms of the behaviour of the mode around the directions of greatest concentration 
of escaping flux, as one might expect intuitively. Further details of this approach are 
developed in Sec. 4. 
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3. Approximation of Herring's integral for generically perturbed cavities 

In this section we describe how the weak-deformation approximation (5) of the Herring 
integral may be exploited to describe the degradation of the quality factor for boundary 
shapes of the form 

r(0) =r + en(d). 

Recall that natural boundaries prevent a systematic calculation of escaping rays in 
generically deformed cavities but that, for small enough deformations, approximate 
ray families calculated using canonical perturbation theory can successfully describe 
the emitted field. We now summarise the main relevant results from that calculation 
and use them to evaluate the integral M. in (5). 

3. 1 . Calculation of external field using perturbative approximation of rays 

The first step in finding the complex action function S'(x) of escaping rays is to 
calculate the corresponding action function on the boundary <9f2. This must be done 
for complex coordinates to allow a continuation to the farfield by tracing the associated 
complex escaping rays. 

The action is found on d£l by applying canonical perturbation theory to the 
internal dynamics, following the approach developed by Prange and Zaitsev in the 
context of billiard problems [24, 25]. Dynamics on the boundary is characterised by 
using a perturbative expansion of the chord function 

L(e,e') = L {e,e') + ^(6, $') + ■■-, (8) 

expressing, as a series in the perturbation parameter e, the length of the chord 
connecting points on the boundary with polar angles 9' and 9. The leading term 

L o (0/)=2r o sin ^ 

is the chord function for a circle of radius r$ and the first-order term is 

'-0' 
2 



Li(0,0') = sin 



(n(6) + n(e')). 



The chord function L(9, 9') serves as a type-one generating function for the boundary 
map, expressed in terms of the polar angle 9 and its conjugate momentum variable 
J = pg. Note that (9, J) serve as action-angle variables for the unperturbed boundary 
dynamics. The action-angle variables (9, J) for the perturbed system can then be 
described by the type-two generating function 

S(0,J) = $J + eg(O,J) + ---, (9) 

where g(9, J) is shown to satisfy the difference equation 

g{0 + u,J)-g{6,J) = R{0,J), (10) 

where 

R(6,J) = L 1 (6 + u,6)-(L 1 ) 

denotes the oscillating part, with respect to 9, of L\(9 + u, 9), with r$ cos(w/2) = J. 
Note that J takes a fixed value determined by the quantisation condition nkJ = m, 
where m is the azimuthal quantum number of the unperturbed mode. 
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This difference equation for g{9, J) is straightforwardly solved, as a Fourier series 
for example, for any analytic function r\{9). Note that the rays escaping to infinity 
typically start with complex values of 9. Here we limit, our perturbative expansion 
to a first-order truncation, in which case it suffices to substitute in g(9, J) the initial 
conditions for the unperturbed escaping rays of the circle. For a counterclockwise- 
rotating mode, these can be shown to originate at6* — x — fi = X ~ k/2 + i0, 
where \ 1S the polar angle of a ray in the farfield and )3 — it/2 — iO is the 
angle of refraction with which it leaves the cavity's boundary. From Snell's law, 
coshO = sin/3 = nsina = nJ/ro, where a is the (constant) angle of reflection of the 
unperturbed internal rays. If the terms in the Fourier series for r\{ff) decay slowly 
enough, the solution g(9, J) of this Fourier series may have natural boundaries which 
lie below the initial conditions on Im(9) = G for escaping rays. In this case it is not 
known how to approximate external field arbitrarily far outside the resonator. We 
therefore restrict our attention to the case where g(9, J) can be successfully extended 
in the complex plane as far as Im(#) = O. This is true in particular when r\{9) is a 
trigonometric polynomial function of 9, which is the case for the models used in our 
numerical illustrations. 

Next, the generating function S(9, J) can also be shown to determine, at first 
order, the phase of the wave field on the boundary. That is, on the boundary, 
ip = B(9)c nlk ^ 0J+e9 ^ B,J ' l+ '"\ The action along a ray emerging from 9 with polar 
angle of propagation \ is then 

5(x, J) = £(x, 9) + nS{9, J) + 0(e 2 ) 

where 

denotes the length of the ray connecting the point with polar angle 9 on dil to the 
point x on S and where, if required, 9 can be expressed in terms of x and J. 

Finally, as the cavity is deformed, the initial coordinate 9 of the ray to a fixed 
exterior point changes. However, the functions £(x, 9) and S(6, J) can be shown 
respectively to satisfy the generating-function conditions dS{9, J)/d9 = r sina and 
dl(x, 6) /d6 = —ro sin/3. By Snell's law, the effect on 5(x, J) of changing therefore 
cancels at first order in e and we may approximate 

5(x, J) = S (x, J) + eS*i(x, J) + 0(e 2 ), 

where 

S 1 (x,J) = t 1 (x,0 o ) + ng(0o,J) (11) 

and <So(x, J) an d #o are ; respectively, the action function and the launching angle for 
the unperturbed limit. 

3.2. Calculating the Q-f actor 

Noting that rays in the farfield of the unperturbed problem satisfy 9o = X~ ft/ 2 + i0, 
the Q-factor degradation (6) can then be approximated 

M = — ic- 2eKlmnx) d X . 
2-k J 

where, 

f(x) = S 1 (x(x),J) 
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Figure 2. The perturbative approximation of k in (12) is compared with a 
numerical evaluation for a deformation of the form r = ro(l + £ cos 30) of the 
circle. The example shown here is a mode with azimuthal quantum number 
m = 60 and 6 radial nodes (with k r ro fb 45.3) in a cavity with refractive index 
n = 2. 



9(0, J) = 



is obtained by evaluating (11) at a position x(x) on T, denned by the polar angle \ 
and using the unperturbed launch angle Oq. 

To illustrate the calculation more explicitly we now consider deformations of the 
form 

n(0) = r cos NO, 
where N is an integer. The difference equation (10) has the solution 

r sin(w/2)sin.ZV6' 
tan(iVw/2) 
in this case, while 

h(x,6) = -r 1 (e)cos( X -0), 
and it can be shown as a result that 

Im/(x) = /ocosiv(x-|) 

where 

nr sinhiVesin(u;/2) . 

/o = 1 — 7T- r /r) x r cosh iV9 sinh 9 

tan(ivw/2) 



Then, 

M=I (2kef ), 
where Iq(z) denotes a modified Bessel function of the first kind. 



(12) 
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A typical implementation of this result is illustrated in Fig. 2. It should be 
emphasised that even though the perturbations involved here are quite small (typically 
of the order of one part in a thousand in a cavity that is about 30 wavelengths in 
diameter), the deformation is enough to significantly change the Q-factor — by a factor 
of 4 or more in the illustration. We also emphasise that, despite the apparent smallness 
of the perturbation, natural boundaries typically prevent us using a straightforward 
exact determination of the rays in this example. 

4. Wilkinson's formula and the ellipse 

For larger deformations, where the relevant eikonal phase is a rapidly-varying function 
of its arguments, Herring's integral (2) invites approximation by the method of steepest 
descents. This approach is attractive because it provides a direct interpretation of the 
Q-factor degradation in terms of simple geometrical properties of particular escaping 
rays — we will see in fact that it takes the form of a formula developed by Wilkinson 
in the context of tunnel splittings of quantum-mechanical energy levels [16]. 

For generic deformations such as treated in the previous section, Wilkinson's 
formula provides us with a simple ray-geometrical interpretation of the large-efc 
asymptotics of approximations such as (12). It should be acknowledged, however, that 
the method will in that case be restricted to moderately large values of ek where the 
perturbative approach remains valid. The asymptotic approach then has the benefit 
primarily of providing a useful ray-geometric interpretation of existing calculations 
rather than pushing them into fundamentally different regimes. 

The case of elliptical deformation, however, is a particularly important special 
case where significantly greater deformations are treatable and where the asymptotic 
approach has the added value of extending the results of the previous section to 
different regimes. 

We begin by describing the appropriate generalisation of Wilkinson's formula to 
the calculation of Q-factor degradation in a more general context where the escaping 
rays (or some approximation of them) are assumed to be known but without making 
further assumptions. We then apply this result to the specific case of elliptical 
cavities where the known invariants of the underlying ray families allow us to provide 
an explicit evaluation of the general result. Finally, we describe at the end of 
this section how a uniform evaluation of Herring's integral allows us to interpolate 
smoothly between the generic perturbative results of the previous section and the 
large-deformation limit described by Wilkinson's formula. 

4-1- Wilkinson's formula in general 

Wilkinson's formula is obtained by carrying out the steepest-descents approximation of 
Herring's integral outlined in section 3. Recall that we assume the existence (exactly, 
as in the ellipse, or more generally following a perturbative approximation) of an 
invariant M of the ray family in addition to the usual ray invariant H = {p\ + pz)/2 
and that the steepest-descents condition (7) is satisfied by a discrete set of real escaping 
rays guiding the maximum-intensity waves to infinity. 

Then expansion of (3) (or, cquivalently, (4)) about each of these stationary points 
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leads to the approximation 

k r d 2 (S-S' 



K = ]T(Rep„)|t| 2 H 



27ri dsds 



1/2 



where 7 labels the real rays escaping to infinity along directions of maximal intensity. 
Further manipulation of the amplitude following techniques given in [26] enables it to 
be expressed in the following canonically invariant form: 

^ ' ' J (14) 






{2nf/ 2 k 1 r /2 y/\{M,M*}\' 



The term in the square root here denotes a Poisson bracket between the invariant 
M and its complex conjugate M*. This Poisson bracket necessarily vanishes on 
real ray families, for which M — M* , so the denominator in Wilkinson's formula 
effectively provides a (canonically invariant) measure of the rate at which rays become 
complex as one moves away from the particular real rays labelled by 7. Alternatively 
it measures the transversality of the ray family to its complex-conjugate partner 
around their real intersection along 7. We also note that, using the fact that M 
and M* each Poisson commute with H = (p 2 + p y )/2, the Jacobi identity yields 
{{M,M*},H} = -{{M*,H},M} - {{H,M},M*} = 0. Because the flow of H is 
along rays, this means that {M.M*} is invariant along each real ray 7 and it therefore 
doesn't matter where on 7 we calculate it. 

4-2. Wilkinson's formula for the ellipse 

As described in [15], the external field can be described in detail in the special case of 
elliptical deformations by exploiting the existence of explicit analytical expressions for 
the invariant M for that case. We now use these results to offer a detailed evaluation 
of Wilkinson's formula for the dielectric ellipse. 

We begin by summarising the main geometric characteristics of the rays of the 
dielectric ellipse, with a more detailed description being available in [15]. We denote by 
a and b the major and minor semiaxes and by c satisfying a 2 + b 2 = c 2 the half-distance 
between foci. Ray families in the ellipse's interior have, as an invariant, 

M*, P) = (Xp y - VPx) 2 - C 2 p 2 y , 

which expresses in cartesian coordinates the product of angular momenta about the 
two foci. The invariant in the exterior takes a different functional form but can be 
expressed in terms of the function A(x, p) by matching interior and exterior families 
using Snell's law. Our convention is that the exterior invariant M thus obtained is 
scaled so that 

n 2 M = A - (n 2 - l)c 2 sin 2 u. (15) 

Here, u denotes the boundary coordinate from which an external ray is launched, where 
in general u parameterises the ellipse boundary according to (x, y) — (ocosu, & sin it). 
Because rays guiding the external wave to infinity have complex starting points on 
the boundary, u and therefore M, are complex functions of the external phase space 
coordinates (x, p). In fact, more detailed calculations [15] show that 

-bp x L + iapyQ 



smti 



b 2 + c 2 p 2 
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where L = xp y — yp x denotes angular momentum about the centre of the ellipse, and 
Q = (A — b 2 (p 2 + Py)) 1 / 2 - We can then separate the real and imaginary parts of the 
invariant M according to 

n 2 M = A- B 2 + C 2 - 2\BC, (16) 

where 

B = -c\Jn 2 - 1 — — x and C = c\/n 2 — 1 — — v . 

b A + c z p~ b A + c z Py 

In order to implement (14) we must now find the real escaping rays 7, along with 
their complex starting points on the boundary, and then evaluate their actions and 
amplitudes in (14). 

Assuming foci on the a;-axis, we find that real rays guiding waves of maximum 
intensity emerge vertically and have starting points with it-coordinates 

ln 2 M -b 2 

u\ = iUi and u 2 = ir + iUi, where sinhC/i = \ -7. tt^t- 

V b z + n z c z 

The full coordinates of these starting points arc (a; ,j/o) — (±acosh?7i, ±ifesinhf7i) 

and the corresponding rays emerge (71 and 72 in Fig. 1) with momentum p = (0, ±1), 

so that a general point on them has coordinates (x,y) — (±acoshUi,yo +t), te C. 

The imaginary part of the action is therefore the imaginary part of the y displacement 

needed to get to real coordinate space, which is 

ImS = Pylm(yo) = bsmhUi. 
It remains to evaluate the Poisson bracket in the denominator of (14), which, from 
(16), takes the form. 

n 4 {M, M*} = Ai{A -B 2 + C 2 , BC}. 

Evaluated on the dominant, vertically escaping rays, this can be shown after further 

manipulation to take the value 

8(n 2 -l)c 2 b, o, ln 2 M-b 2 , N 

i{M,M*} UuU2 = n J (M + c 2 ) \j b2 + n2c2 - (17) 

The final ingredient needed to evaluate (14) is the transmission coefficient, which, 
following the discussion in [15], can be shown to take the form 

2i/n sin a cos B ■ 

t = e 

cos a + cos /? 

where a and /3 are respectively the angles of incidence and reflection at the boundary 

point from which the escaping ray is launched. Our convention in this paper is that 

this transmission coefficient also incorporates a phase correction a, described in detail 

in [15], that accounts for the variable reflection phases on the boundary and is related 

to the Goos-Hanchen shift [22, 23]. 

A numerical illustration of this result is given in Fig. 3. Also shown in the same 

figure is the perturbative prediction developed in the previous section. Approximation 

(14) is seen to give a good description of the quality factor degradation, and it 

works well for deformations strong enough that the perturbative result (12) has 

failed. Approximation (14) fails, however, for very small perturbations where, in the 

intcgrable circular limit in which all escaping rays become real, the Poisson bracket 

{M,M*} vanishes and the amplitude of (14) diverges. This is because in the regime 

of small ke the steepest descents approximation used to approximate n is no longer 

justified. Here the perturbative approach of the previous section must be used instead. 

In the next subsection we describe how the two limits may be married in a single, 

uniform approximation that is valid in both limits. 
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Figure 3. A comparison is given between calculations of k using Wilkinson's 
formula (14), the perturbative approach of the previous section, and a numerical 
evaluation. The mode is the same as used in Fig. 1 and the perturbation parameter 
e is defined so that b = ro(l — e) (so that c 2 = a 2 — b 2 rj ier^). 



4-3. Uniform approximation 

A uniform approximation of the Q-factor degradation that interpolates smoothly 
between the perturbative approximation of Sec. 3 and Wilkinson's formula in (14) 
is now derived. The method here follows a uniform treatment in [27, 28] of the 
contribution of bifurcating periodic orbits to the trace formula. 

We return to the version of Herring's integral given in (4) and seek a coordinate 
transformation s — > \ which is so that the imaginary action takes the form 

ImS = K + AK cos2 X , 

where K and AK are constants. The form chosen here exploits the symmetry of the 
problem with respect to u — > u + 7r and the variable \ thus defined coincides at leading 
order with the polar angle used in the perturbative Herring integral (5), but deviates 
from it at higher deformation. 

The constants K and AK are determined by evaluating the minimum and 
maximum values K m [ n = K — AK and K max = K + AK of ImS*. We have already 
established that the minimum value K m - m = b sinh U\ is achieved by the real rays 
escaping along the directions of greatest intensity. The maximum value if ma x achieved 
by real rays (73 and 74 in Fig. 1) escaping horizontally along directions of least 
intensity. It can be shown that these latter rays are launched from points 



u = ± — h 1U2 , where cosh XJ<z = 
and have an imaginary action Im S — K max = a sinh [7 2 ■ 



n 2 M 
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Following this change of coordinates, the Herring integral is left in the form 

K = — (^ A(x)e- k " ( - R+AKcos2 ^d X . (18) 



We next write 

Mx) = A v + A x cos 2 X + H( X )g-AK cos 2 X , (19) 

where Aq and A\ are constants, chosen so that 

ff( . Ap + Ai- A{ x ) 
{X> 2AKsin2x 

is smooth, which (exploiting the symmetry of the problem with respect to u — > u + it) 
is true provided 

1(0) = A Q + A x (20) 

and 

A(*)=A -A 1 . (21) 

Substitution of (19) in (18), allows us, following integration by parts, to neglect at 
leading order in l/k r the contribution of the last term in (19). The result is 

k w (A I (k r AK) + A 1 I 1 (k r AK))c- k '' R , (22) 

where Iq{z) and I\(z) denote modified Bessel functions of the first kind. 

In practice, the amplitudes Aq and A\ are conveniently evaluated by comparing 
the large-deformation asymptotics of (22) with the amplitude terms of Wilkinson's 
formula. A straightforward comparison yields 

{A , A) H^Z = l l^)l 2 ^ 

1 0+ '^KAK ( 2T )3/2 fc J/V|{M,M*} tll | 

Although the physical problem insists that k r is positive, we might also formally 
compare asymptotic results for negative k r and this yields 

1 \t(u 3 )\ 2 J 



(A - A,) 



k r AK (2^)3/2^/2 ^/\{M,M*} U3 \ ' 
The right-hand side here is determined by the geometry of the ray family around real 
rays escaping along the directions of least emitted intensity. Here we need the dark 
analogue of (17): 



. rw „«-, 8(n 2 -l)c 2 a r n 2 M -a 2 +n 2 c 2 ,„„, 

n A b z V a — n A c A 

Aq and A\ are then known and the evaluation of (22) is complete. 

This result is compared with numerical data in Fig. 4 and is seen to give a very 
good description of the Q-factor degradation from the perturbative regime of (12) 
to the larger-deformation regime of (14). In practice, the deformations treatable by 
this result are limited by our ability (or otherwise) to find phases a — which have 
been absorbed in the transmission factor \t\ 2 in this paper — which account for the 
dielectric boundary conditions in this problem. Such phases have been seen in [15] to 
suffer from natural boundaries which approach the real axis as eccentricity increases 
and invalidate the WKB ansatz above a critical deformation in the region of the 
complex boundary from which escaping rays are launched. It should be emphasised 
that the eccentricities at which this happens are formally of 0(1) and much greater 
than the eccentricities at which the perturbative result fails. 
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Figure 4. A comparison is given between calculations of n using the uniform 
result (22) and a numerical evaluation for the same mode and parameter ranges 
used in Fig. 3. 



5. Conclusions 



We have found that, in a regime of short wavelength, very small deformations of optical 
cavities suffice to alter significantly the quality factor of resonant modes and have given 
quantitative estimates of this effect. We have shown that the quality factor may be 
halved or worse by deformations as small as one part in a thousand, for example, in 
explicit calculations where the cavity is some tens of wavelengths in diameter. 

Although complex WKB methods provide a natural starting point for such 
calculations, their implementation is problematic because the required ray data may 
be difficult to find, or may not even exist. Instead, we have provided a general analysis 
based on approximation of the underlying ray families using perturbation techniques. 
Although restricted to relatively weak deformations, because the perturbations are 
applied to the rays rather then directly to the wave solution itself, this approach can 
nevertheless successfully describe deformations where the wave and quality factor itself 
are altered nonperturbatively. 

For moderately large deformations, the results of this analysis are naturally 
interpreted in terms of the geometry of families of escaping rays, using Wilkinson's 
formula (14). Deformation effects the important qualitative change of making the 
associated family of escaping rays slightly complex (whereas far enough away from the 
cavity the corresponding orbit families are real in the undeformed case). Nevertheless 
a discrete subset of the escaping rays remain real and these in particular determine 
the directions along which the waves of greatest intensity propagate. Wilkinson's 
formula expresses the quality factor degradation in terms of the geometry of the ray 
family around the real subset guiding waves of highest intensity. It fails in the limit 
of small deformation where a denominator, which measures the rate at which rays 
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become complex away from the real subset, vanishes. A uniform result has been 
also derived, however, which interpolates between the primitive Wilkinson formula 
appropriate to larger deformations and the perturbative analysis derived earlier for 
weak deformations. This was shown to give an excellent description of the quality 
factor in the special case of elliptical deformations, where the ray families can be 
calculated analytically even for significantly eccentric cavities. 
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